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Abstract —In this paper, we study the average throughput 
performance of energy beamforming in a multi-antenna wireless- 
powered communication network (WPCN). The considered net¬ 
work consists of one hybrid access-point (AP) with multiple 
antennas and a single antenna user. The user does not have 
constant power supply and thus needs to harvest energy from 
the signals broadcast by the AP in the downlink (DL), before 
sending its data back to the AP with the harvested energy 
in the uplink (UL). We derive closed-form expressions of the 
average throughput and their asymptotic expressions at high SNR 
for both delay-limited and delay-tolerant transmission modes. 
The optimal DL energy harvesting time, which maximizes the 
system throughput, is then obtained for high SNR. All analytical 
expressions are validated by numerical simulations. The impact of 
various parameters, such as the AP transmit power, the energy 
harvesting time, and the number of antennas, on the system 
throughput is also investigated. 

Index Terms —RF Energy harvesting, energy beamforming, 
wireless-powered communication network, performance analysis, 
throughput. 


I. Introduction 

Limited energy supply has been a serious problem faced 
by many wireless communications systems due to a dramatic 
increase of their energy consumption. Recharging and replac¬ 
ing batteries may introduce a high cost and sometimes could 
be inconvenient or impossible, as for example, in dangerous 
and toxic environments or for sensors underwater or inside 
the human body. In this context, energy harvesting techniques, 
which can prolong the lifetime of energy constrained wireless 
networks in a sustainable way, have emerged and attracted 
significant attentions, e.g., see H] and references therein. The 
techniques enable wireless nodes to harvest energy from the 
surrounding environment. Apart from other techniques that 
harvests energy from solar, wind, vibration, thermoelectric 
effects or other physical phenomena, radio frequency (RF) 
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energy harvesting, which typically refers to the capability 
of the wireless nodes to scavenge energy from RF signals, 
has emerged as an attractive solution recently El, El- The 
feasibility of RF energy harvesting has been experimentally 
demonstrated by hardware implementation in HI and 0. 
Wireless-powered communication networks (WPCNs), a new 
type of wireless networks emerging along with the RF energy 
harvesting technique, refer to networks in which the wireless 
devices normally have no internal energy supply (e.g., battery) 
and need to harvest energy from the RF signals emitted 
from a (dedicated) power transmitter before their information 
transmission Jb]. In practice, WPCNs can find many potential 
applications, such as wireless sensor and RFID networks Q. 

There have been several research papers in open literature, 
which studied WPCNs for different network setups. Single- 
user WPCNs were investigated in 0, El- Taking into con¬ 
sideration limited feedback, Chen et al. derived the optimal 
time duration for the downlink (DL) wireless energy transfer 
(WET) through maximizing upper and lower bounds on the 
average information transmission rate of one-user WPCNs 
0. Furthermore, the impact of channel estimation error on 
the average information transmission rate was also analyzed. 
0 optimized the energy efficiency of a single-user WPCN 
by jointly designing the time duration and transmit power 
for DL wireless energy transfer. Multi-user WPCNs were 
first investigated in GO), in which a “harvest-then-transmit” 
protocol was developed. In the proposed protocol, the users 
first collect energy from the signals broadcast by a single¬ 
antenna hybrid access-point (AP) in the DL and then use 
the harvested energy to send independent information to the 
hybrid AP in the uplink (UL), based on the time-division- 
multiple-access (TDMA) scheme. The sum-throughput of the 
considered network, subject to user fairness, was maximized 
through optimizing the time allocation for DL WET and UL 
wireless information transfer (WIT). A similar network setup 
with a multi-antenna AP was investigated in ED, where 
multiple users can simultaneously transmit information to the 
AP in the UL through space-division-multiple-access after they 
harvest energy in the DL. In ATI , the minimum throughput 
of all users was maximized by jointly optimizing the time 
allocation, the DL energy beams, the UL transmit power 
allocation, and the receive beamforming vectors. 

In contrast to the above papers, which mainly focused on 
resource allocation in WPCNs, in this paper we concentrate on 
analyzing the average performance of one-user multi-antenna 
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WPCNs, in order to gain important insights for designing and 
implementing multi-antenna WPCNs in practice. Specifically, 
we derive the exact closed-form expression for the average 
throughput performance of a multi-antenna WPCN with en¬ 
ergy beamforming. It is also worth mentioning that there are 
several papers that analyzed the average performance of the si¬ 
multaneous wireless information and power transfer (SWIPT) 
for different network scenarios mi-iiisi, which are related 
but are essentially different to this paper. These papers focused 
on the characterization of the fundamental trade-offs between 
energy transfer and information transmission using the same 
signal. Moreover, im-iiiii focused on the development of 
cooperative protocols for WPCNs with different setups. Note 
that all nodes were assumed to be equipped with single antenna 
in these papers. To the authors’ best knowledge, this is the 
first paper that derives exact average performance of a multi¬ 
antenna WPCN with energy beamforming. 

In this paper, two different transmission modes are consid¬ 
ered for the user, i.e., delay-limited and delay-toleranH which 
corresponds to different length of the code-words used by the 
user in. Col. More specifically, in the delay-limited mode, 
the received signal at the AP has to be decoded block by block 
and the code length should be no longer than the time duration 
of each transmission block. Thus, the average throughput of 
this transmission mode can be obtained via evaluating the 
outage probability. In contrast, in the delay-tolerant mode, the 
AP can store the received information blocks in a buffer and 
tolerate the delay for decoding the stored signals together. In 
this case, the code length can be very long compared to the 
transmission block time. The egordic capacity is then used to 
calculate the system average throughput in this delay-tolerant 
mode. 

The main contributions of this paper are summarized as 
follows: 

• We derive exact closed-form expressions for the average 
throughput of the considered network in both delay- 
limited and delay-tolerant transmission modes by eval¬ 
uating the outage probability and ergodic capacity with a 
given energy harvesting time. 

• Asymptotic analyses are performed to gain insights and 
derive closed-form expressions of the throughput-optimal 
energy harvesting time at high SNR for both transmission 
modes. Futhermore, numerical simulations are presented 
to provide practical design insights into the impact of 
various parameters, such as the transmit power of the AP, 
the energy harvesting time, and the number of antennas, 
on the system performance. 

Notations: ||x|| denotes the Euclidean norm of a vector 
X. Fx (x) and fx (x) are used to represent the cumulative 
distribution function (CDF) and probability density function 
(PDF) of a random variable X, respectively. 

The rest of the paper is organized as follows. Section 
II describes the system model. The exact and asymptotic 
performance analysis of the considered system is provided in 
Section III. Section IV presents the numerical results from 

^The delay-tolerant mode is named as no-delay-limited in (20). 
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Fig. 1. System model of a one-user multi-antenna WPCN. 

which various design insights are given. Finally, Section V 
concludes this paper. 

II. System Model 

We study a single-user multi-antenna WPCN shown in Fig. 
[T] with energy transfer in the downlink (DF) and information 
transmission in the uplink (UF). We follow 0, mil and 
consider that the AP is equipped with N antennas but the 
user carries one single antennqS Such a network setup could 
correspond to a wireless sensor network in practice, where the 
user is a sensor node that cannot be equipped with multiple 
antennas due to the size and cost constraints. Both the AP and 
user work in a half-duplex mode. It is assumed that the user 
has no embedded energy supply and thus needs to first harvest 
energy from the signals broadcast by the AP in the DF and 
then use it for the information transmission to the AP in the 
UF. In our model, we consider a battery-free energy storage 
design at the user, which has been adopted in several recent 
works such as 0, ifTOll . ifT^ . ifTTll . Specifically, the user just 
uses a storage device such as a supercapacitor to keep the 
harvested energy among each transmission block. 

The “harvest-then-transmit” protocol proposed in mni is 
assumed to implement in the considered network. Fet T denote 
the duration of one transmission block. The first interval of 
the duration tT in each transmission block is used for the DF 
wireless energy transfer from the AP to the user and in the 
remaining interval of the duration (1 — t)T, the user applies 
the harvested energy to transmit the information to the AP 
in the UF . For simplicity but without loss of generality, we 
consider a normalized unit block time (i.e., T = 1) hereafter. 
In line with 0, ifTTl . we assume a quasi-static channel model 
with perfect channel state information (CSI) at the AlH In 
practice, the CSI can be acquired by various methods, e.g., 
the pilot-assisted reverse-link channel training I^ . It is worth 
mentioning that it would be more general to consider the 

^The extension to a general m ultip le -inp ut multiple-output (MIMO) sce¬ 
nario with single/multiple user(s) CS, Eil would require the development 
of a new framework to analyze its system performance, which may constitute 
an independent paper and thus out of the scope of this paper. 

^In this paper, we consider the case that at the beginning of each trans¬ 
mission block, the AP first transfers a certain amount of energy to the user 
to acquire CSI. After this CSI acquisition duration, the AP delivers wireless 
energy to the user again, which will be used for UL information transmission. 
We assume that the channel is either static or varies very slowly, which is the 
case for wireless sensor networks, and thus the CSI acquisition duration is 
very short such that it can be ignored compared with the whole transmission 
block. 
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imperfect channel model as in i). But this would require 
a totally different performance analysis framework, which 
may constitute an independent paper. Thus, we would like 
to consider this as our future work. 

Let h and g represent the N dimensional DL and LfL 
channel vectors, respectively. Their entries are assumed to 
follow independent and identically distributed (i.i.d) circularly 
symmetric complex Gaussian distribution^ with a zero mean 
and variance of fl. Let P denote the transmit power from 
the AP during the DL phase. With the CSI h available at the 
AP, the maximum ratio transmission (MRT) with beamforming 
vector w* = h/||h|| is known to be optimal for DL energy 
transfer il. By ignoring the negligible energy harvested from 
the noise, the amount of energy harvested by the terminal 
during the period r, denoted by Eh, is given by 

Eh=TlTP\\H\ ( 1 ) 

where p is the energy conversion efficiency at the user. 

After harvesting the energy from the AP, the user transmits 
the information to the AP during the remaining time 1 — r. 
The received signal at the AP can be expressed as 


yA = 



gs + n, 


( 2 ) 


where s is the transmitted signal, and n is the additive Gaus¬ 
sian white noise with a zero mean and variance matrix o^In- 
The maximum ratio combing (MRC) technique is assumed to 
be performed at the AP to maximize the received signal-to- 
noise ratio (SNR). The resulting SNR at the AP can be written 
as ||9] 


EhM? 


(1 — r) (T^ 



(3) 


where 7 = 


III. Performance Analysis 

In this section, we analyze the average throughput perfor¬ 
mance of the considered system in both delay-limited and 
delay-tolerant transmission modes, respectively. This will be 
achieved by evaluating the outage probability and ergodic 
capacity of this network. We will also perform asymptotic 
analysis to obtain some insights on the effect of various 
system parameters. Moreover, we derive the expressions of 
the throughput-optimal energy harvesting time r at high SNR 
for both transmission modes. 


A. Delay-Limited Transmission Mode 

1) Exact Throughput Analysis: In this transmission mode, 
the throughput can be obtained by evaluating the outage 
probability defined as the instantaneous channel capacity drops 
below the source’s fixed transmission rate R IfT^ . Thus, the 


^Note that Rice fading is another practical model to describe the channel 
fading between the AP and user when the distance between them is limited and 
the light of sight exists. The performance analysis of the considered network 
over Rice fading has been treated as our future work. 


outage probability of the considered system with a given rate 
R and energy harvesting time r can be expressed as 


Poutir) = Pr (log 2 (1 -b 7 a) < R) = Pt ( 7 a < 7o) 

/ I -\ 

(\/W) 


= 1 - 


N-l 

2 E 

p—0 


Kn-p 2 


T)7'n 1 


(iV-l)! 


(4) 


where 70 = 2^ — 1 and Kn (•) is the modified Bessel function 
of the second kind with order n ll25l Eq. (9.6.2)]. The last 
equation is obtained by noting that the terms ||h||^ and ||g||^ in 
7 A follow independent and identical chi-square distributions 
with 2N degree of freedom ED, and applying the PDF of 
product of two chi-square distributed variables derived in Il 22 l . 

Considering that the user transmits with a fixed rate R and 
the effective communication time from the user to AP after 
the energy harvesting phase is 1 — r, the average throughout 
of the delay-limited transmission mode, pum, is given by 


Plim = (1 - Poutir)) R{1 -t) 


2R{1 - t) 


N-l 

E 

p—O 


_\ N-\-p 

(l-r)70 \ 






N-p 


9 / ( 1 - T ~)70 


p\ 


(5) 


where the throughput pum depends on r, N, 70 , Pi, and 7 . 
Although the expression Q gives exact value of the system 
throughput and may not be computationally complicated to 
evaluate, it does not offer explicit insights into the impacts of 
the system parameters. Thus, we are motivated to look into the 
system performance for this transmission mode at high SNR, 
where simpler expressions can be obtained. 

2) Asymptotic Throughput Analysis: To this end, we can 
apply the series representation of Bessel functions in (|5]l. But, 
due to complex structure of (jS), this method is shown to be not 
tractable mathematically. To address this, let us first derive a 
tight approximation for (|5]). Recall that the end-to-end SNR 7 a 
in (O can be regarded as a product of two chi-square random 
variables. In addition, the chi-square distribution can be re¬ 
garded as a special case of the generalized gamma distribution. 
Recently, Chen et al. proposed a novel analytical framework 
for evaluating the statistics of the product of independent 
random variables in EH, where a tight approximation to the 
CDF of the product of generalized gamma random variables 
was provided. By applying the results obtained in Ell, we 
can approximate the CDF of 7 a as 


PlA (^) 


7(mo + 2iV-2,^V^) 
r (mo + 2N -2) 


(6) 


where y{a,x) = is the lower incomplete 

gamma function and F (a) = is the complete 

gamma function. Besides, mg = 1.6467 and fig = 1.5709 are 
heuristically obtained ll24l . 

Consequently, the average throughput of the delay-limited 
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transmission mode can be approximately expressed as 


Pli 






r (mo + 2iV - 2) 


(7) 

As shown by numerical results in Sec. IIVI the above approx¬ 
imation is very tight. 

In high SNR regime (i.e., 7 —>• 00 ), the term 
inside the incomplete game function approaches to 0. By 
applying the asymptotic property of the incomplete gamma 
function near zero 7 (a, x) k, x°'/ a, we further have 


Pli 




1 - 


2mQ / (1 -t)7q \ 

Oq Y j 


_\ mo+2A^ —2' 


r (mo + 2N — 1) 


( 8 ) 


From we now can observe the impact of the system 
parameters on the throughput performance. Specifically, when 
any of the parameters N, Q and 7 increases, the ratio term in¬ 
side the square brackets o/® will decrease, thereby improving 
the system throughput. However, the system throughput cannot 
keep increasing because the value of the term inside the square 
brackets is limited by 1 no matter how large are the values 
of N, n and 7 . In other words, the system throughput in the 
delay-limited transmission mode should converge to a ceiling 
value of i?(l — t) by increasing the parameters N, and 
7 . Moreover, the impact of the parameter r is not so explicit 
as that of the aforementioned ones. Particularly, the value of 
the term inside the square brackets rises when the value of 
T increases. But, this will lead to the decrease of the term 
(1 — r) at the same time. Thus, we claim that there should 
exist an optimal energy harvesting time, denoted by that 
can achieve the maximal throughput of the delay-limited mode. 

Based on the approximate system throughput ®, we can 
derive an approximate expression for the value of in high 
SNR regime, which is given in the following proposition: 

Proposition 1: The throughput-optimal energy harvesting 
time at high SNR in the delay-limited transmission mode can 
be approximately expressed as 



B. Delay-Polerant Transmission Mode 


1) Exact Throughput Analysis: In this subsection, we an¬ 
alyze the average throughput performance of the considered 
system with a delay-tolerant transmission mode. In this trans¬ 
mission mode, it is assumed that the AP can store the received 
information blocks in a buffer and tolerate the delay for decod¬ 
ing the stored signals together. Then, the average throughput 
can be determined by evaluating the ergodic capacity at the 

AP mi. 

We first evaluate the ergodic capacity of the considered 
network for a given energy harvesting time t. According to 
the definition, the ergodic capacity of the system, denoted by 
C (t), is given by 


pCO 

C (r) = / log 2 (1 + 2 ) {z)dz 
Jo 


y J 


N 


(iV- l)!ln2 (iV- 1)! 

^ \n{l + z)z^-^Ko(2^ 


( 10 ) 




dz. 


where (z) denotes the PDF of the random variable ja, 
which can be expressed as 


flA (Z) = 




77 F j 


,Ar-l 


Kn 


n-r)z 


( 11 ) 


(iV- l)!(iV- 1)! 


Note that the ergodic capacity can also be calculated by using 
the CDF, which is, however, essentially the same with the 
adopted method using PDF. 

Inspired by ll2^ . we can evaluate the integral term in (flOl) in 
a closed-form by applying the Meijer G-functions as follows. 
Firstly, we express the terms In (1 -f a;) and x^'K^, (x) in 
Meijer G-functions by using the equations given in lIZTl Ch. 
2]. Then, we can solve the integral and write (fTOl) in a closed- 
form through the following calculation 


C{r) = K^x 




rGiiUz 


1,1 

1,0 


G: 


!,0 f {1 -t)z 

),2 Y 


N -l,N -l]dx 


l — T ,^4,1 j l — T 

r^n^^2,4 I 'rfrP 


(iV- 1)!(W- l)!ln2 

- 1,0 

-1,N - 1 


(Af- 1)!(A^- l)!ln2 


( 12 ) 


where W(x) is the Lambert W function, which is the solution 
of the equation W e'^ = x. Besides, A = mo -f 2N — 2 and 
B = \J defined for notation simplification. 

Proof: See Appendix lAl ■ 

Note that the above expression for the optimal energy 
harvesting time is derived for high SNR regime. To obtain 
the optimal r of any given SNR, we can calculate the first- 
order derivative of ® with respect to r and set it equal to zero. 
However, due to the complexity of ®, it is hard to achieve 
closed-form expressions for the roots of the obtained equation, 
which can only be calculated via numerical methods. 


where G™’" 


is the Meijer G-function 


Qj\ , . . . , Qp 
bi,...,bg 

Eq. (9.301)], which is a standard built-in function in most of 
the available mathematical software packages, such as MAT- 
LAB, MAPLE and MATHEMATICA. Besides, the integral 
of Meijer G-function is solved based on ll^ Eq. (7.811)]. 
Now, we can obtain the average throughput of the considered 
system, denoted by ptoi, by calculating the product of the 
ergodic capacity G(r) and the effective time duration for 
information transmission of 1 — r. That is. 


Ptoi = (1 - r) G (r), 


(13) 
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where the average throughput ptoi depends on r, N, fl and 
7. However, due to the complexity of the Meijer G-functions, 
it is hard to observe the specific relationships between the 
throughput and the aforementioned parameters from (fTSl l. 

2 ) Asymptotic Throughput Analysis: To gain more insight, 
next let us derive the asymptotic throughput. For sufficiently 
high SNR, we can approximate log 2 (1 + z) in the integral of 
(fTOl i as log 2 (z). The asymptotic throughput can be calculated 
without the need of Meijer G-function as 


Ptol 


f - \^ 

2(l-r) \T^) 

{N -ly. (iV- l)!ln2' 


z^-^nzKo 2 ^ 


2(1 — t) j 

{N -ly (iV- l)!ln2 

^ Ko(2 

1 — T 


' (1 — r) z 


dz 


\/P') 


In 2 


2'!/; (iV) + In 7 + 2 In H — In 


2tdt 
1 — T 


(14) 


where Il28l Eq. (2.16.20.1)] is used to solve the integral, ip (•) 
is the Euler Psi function ll23l Eq. (8.36)]. From ( I74l) . we can 
see that for a given value of t, the system throughput is 
proportional to the Psi function of the number of antennas 
N at high SNR. This indicates that the system throughput 
increases with the rising of the number of antennas but 
the increasing rate is gradually decreasing, according to the 
properties of the Psi function. The average throughput is also 
proportional to the logarithm functions of the parameters 7 
and H in a high SNR regime. In addition, similar to the case 
in the delay-limited mode, we can clearly observe from m 
that the energy harvesting time r plays two opposite roles in 
the system throughput. Specifically, increasing the value of r 
will decrease the value of the term outside the square brackets 
but increase the value of the term inside the square brackets. 
Thus, we also deduce that there exists an optimal value for 
the energy harx’esting time, denoted by that can maximize 
the system throughput. 

Based on (0, we can obtain a closed-form expression 
for the optimal energy harvesting time at high SNR for 
the delay-tolerant transmission mode, given in the following 
proposition: 

Proposition 2: The throughput-optimal energy harvesting 
time at high SNR in the delay-tolerant transmission mode can 
be expressed as 


1 -f W ■ 

Proof: See Appendix iBl ■ 

It is worth emphasizing that the above equation for the 
throughput-optimal energy harvesting time is actually a high- 
SNR approximate solution as it is obtained based on the high- 
SNR approximation in (fT4l) . 


From (O, we can see that the value of the optimal energy 
harvesting time t* is inversely proportional to the parameters 
7, n and N, since the Lambert W function W(x) is a 
monotonically increasing function for x > 0. This observation 
will be verified by simulation results in next section. 

IV. Simulation Results 

In this section, we provide simulation results to validate the 
derived analytical results. We also investigate the effects of 
the transmit power of the AP P, the number of antennas N, 
and the energy harvesting time r, on the system throughput. To 
capture the effect of path-loss on the network performance, we 
use the channel model that H = 10~^ {dAu)~°‘, where dAu 
is the distance between the AP and the user, and a € [2,5] 
is the path-loss factor 1^ . Note that a 30dB average signal 
power attenuation is assumed at a reference distance of Im 
in the above channel model nni. In all following simulations, 
we set the distance between the AP and user dAu = 10m, the 
path-loss exponent a = 2, the noise power = —80dBm, 
the energy harvesting efficiency rj = 0.5, and the transmission 
rate of the user i? = 2 in the delay-limited mode. 

Pig.|2]plots the average throughput curves of the considered 
network versus the transmit power of the AP for both delay- 
limited and delay-tolerant modes with various numbers of 
antennas at the AP. We can see from Pig. |2] that the derived 
exact analytical results tightly coincide with the corresponding 
simulated ones, which verified the correctness of the derived 
closed-form expressions for the average throughput given in 
© and (Ql. Besides, it can be observed from Pig. |2](a) 
that the approximate throughput of the delay-limited mode 
calculated via © can estimate the exact one very well, 
which confirms the effectiveness of the approximation. More¬ 
over, for the delay-tolerant mode, the asymptotic throughput 
quickly approaches the exact one as the AP transmit power 
increases in both transmission modes, although there are 
gaps between them when the transmit power is low. This 
observation validates our asymptotic analysis in (fT4l i. It can 
also be observed from this figure that the system throughput 
in both transmission modes is improved as either the transmit 
power or the number of antennas increases. However, the 
trends for the increasing of the system throughput are different 
in two transmission modes. Specifically, in the delay-limited 
transmission mode, the system throughput tends to be saturated 
when the transmit power of the AP is high enough. In addition, 
the larger the number of antennas at the AP, the lower the value 
of the transmit power from which the saturation of the system 
throughput starts. In contrast, the system throughput grows 
monotonically as the transmit power of the AP increases. 
These observations are actually consistent with our theoretical 
analyses discussed in Sec. |III1 

Pig. [3 illustrates the effect of the number of antennas at the 
AP on the system throughput, in which the throughput curves 
of both transmission modes are plotted versus the number of 
antennas for different values of the energy harvesting time 
r. As we expected, the system throughput in delay-limited 
mode first increases and then converges to a ceiling as the 
number of antennas at the AP increases. This is because that 
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(a) Delay-limited mode (b) Delay-tolerant mode 

Fig. 2. Average throughput versus the transmit power of the AP for (a) delay-limited transmission mode and (b) delay-tolerant transmission mode with 
different number of antennas, where r = 0.5. 




(a) Delay-limited mode (b) Delay-tolerant mode 

Fig. 4. Average throughput versus the energy harvesting time r for (a) delay-limited transmission mode and (b) delay-tolerant transmission mode with 
different values of N and P. 



Fig. 3. Average throughput versus the number of antennas at the AP for 
both delay-limited transmission mode and delay-tolerant transmission mode 
with different values of r, where P = 30 (dBm). 


the system throughput in this delay-limited mode is limited 
by the value of i?(l — r). On the other hand, the system 
throughput increases accordingly when the number of antennas 
becomes larger in the delay-tolerant mode. However, the slope 
of the throughput curves gradually decreases, which coincides 
with our theoretical analysis because the average throughput 
is proportional to the Psi function of the number of antennas. 
Moreover, we can see from Fig. [3] that the impact of the 
energy harvesting time r on the system throughput of both 
transmission modes is not so clearly as that of the number of 
antennas does. For example, the throughput performance in 
the delay-tolerant mode is improved when the value of r is 
increased to 0.3 from 0.1. However, as t increases from 0.3 to 
0.5, the system throughput for N >2 significantly decreases. 

To clearly demonstrate the impact of the energy harvesting 
time r, we plot the throughput curves versus r for both trans¬ 
mission modes with variable numbers of antennas and transmit 
powers of the AP in Fig. |4] It can be observed from this 
figure that there exists a throughput-optimal energy harvesting 
time for all the considered cases in both transmission modes, 
which coincides with our analytical results. Recall that we 
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also derive the closed-form expressions for the optimal energy 
harvesting time at high SNR for both transmission modes in 
Sec. ini] To validate these two derived expressions given in (|9]l 
and (ffSl) . we compare them with the corresponding optimal 
energy harvesting time found via exhaustive search in Fig. 
|5] We can see from Fig. |5] that the values of the optimal t, 
calculated by (|9l) and ( fT5l l. coincide with the corresponding 
ones obtained via exhaustive search when the transmit power 
of the AP is high enough (i.e., at high SNRs), which verifies 
the correctness of our derivations of these two equations. In 
addition, as expected, the values of the optimal r decrease as 
either the number of antennas or the transmit power increases. 
This is because the user can harvest the same amount of energy 
in a shorter time when either the number of antennas or the 
transmit power increases, and more time could be allocated 
to the UL information transmission to improve the system 
throughput. 


V. Conclusion 

In this paper, we derived exact and asymptotic closed- 
form expressions for the average throughput of a multiple- 
antenna wireless-powered communication network with en¬ 
ergy beamforming in both delay-limited and delay-tolerant 
transmission modes. The correctness and effectiveness of 
these theoretical analysis were verified by simulation results. 
Numerical results showed that the system throughput in the 
delay-limited mode first improves as either the number of 
antennas or the transmit power of the AP increases and then 
get saturated when these two parameters are large enough. 
In contrast, the system throughput in the delay-tolerant mode 
is always proportional to the number of antennas and the 
transmit power of the AP. Based on the asymptotic analysis, 
we also derived closed-form expressions of the throughput- 
optimal energy harvesting time at high signal-to-noise ratio 
(SNRs) for both transmission modes. Numerical simulations 
showed that the value of the optimal energy harvesting time 
is inversely proportional to the number of antennas and the 
transmit power of the AP. The optimal energy harvesting time 
calculated by the derived expressions can coincide well with 
the one obtained via exhaustive search when the SNR is high 
enough. Furthermore, the larger the number of antennas, the 
more accurate the optimal energy harvesting time calculated 
by the derived expressions. 

Appendix 

A. Proof of Proposition Q] 

The approximate expression of the system throughput given 
in 0 instead of the exact one is adopted to calculate the 
optimal energy harvesting time since it is more tractabl^ 
To proceed, we calculate the first-order derivative of 0 with 
respect to r and set it equal to 0. After some algebraic 

^Originally, the asymptotic expression should be used to find the optimal 
energy harvesting time. However, due to the special polynomial structure of 
we found it is hard to obtain a closed-form expression for the optimal 
energy harvesting time from i). Thus, we choose to use 0 instead. 


manipulations, we have 


A,B 


1 — T 


1 


B 


1 — T 




= r(A), 


and B = 2^ 



where A = mo + 2N — 2 
high SNR, the term B approaches to 0. By applying the 
asymptotic property of the incomplete gamma function near 
zero 7 (a,x) ~ x°‘/a once again and performing variable 


substitution y = 


can simplify (fTbl l as 




(17) 


Furthermore, in high SNR regime, only a small fraction of 
each transmission block is needed for DL energy transfer and 
the majority of time should be allocated to the UL information 
transmission. This means that y = \J~^^ ^ 1 should hold. 
In this case, we can ignore the two “1” terms on the left-hand 
side of (fTTl i and obtain the following equation 


A+2 -By _ (^) 

^ ~ B^ ■ 


(18) 


By equaling the {A-|-2}th roots of both sides in (fTSl l. we have 

'2r(A)' 


ye 
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■ A+2y — 
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A + 2 
B 

A + 2 


ye 


B^ 


T^y = 


B f2r{A) 


y = W\- 


A -f 2 V 

B f2T{A) 


(19) 


A + 2\ B^ 


Substituting the relationship y = into the last equation 

and rearranging it, we can obtain the desired results given in 

0. 


B. Proof of Proposition \2\ 

To find the optimal energy harvesting time t*^j, we calculate 
the first-order derivative of the right-hand side of (fT4l) with 
respect to r and set it to zero. After some algebraic manipu¬ 
lations, we obtain the following equation: 

2^1 (N) + In --=\n . (20) 

T T 

By performing an exponential operation on both sides of the 
above equation, we have 

T 

r 

= ( 21 ) 

r 

= w , 

1 

1 + w ’ 

which completes the proof. 


























(a) Delay-limited mode (b) Delay-tolerant mode 

Fig. 5. The optimal energy harvesting time versus the transmit power of the AP for (a) delay-limited transmission mode and (b) delay-tolerant transmission 
mode with different values of N. 
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